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Abstract 

We study the gravity assist in the general case, i.e. when the spacecraft is not in a coplanar 
motion with respect to the planet's orbit. Our derivation is based on Kepler's planetary motion 
and Galilean addition of velocities, subjects covered in introductory physics courses. The main 
purpose of this paper is to illustrate how the gravity assist can be used to deviate a spacecraft 
outside its original plane of motion. As an example, we use the NASA-ESA's Ulysses mission to 
"test" our simple model. 
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I. INTRODUCTION 



It is well known that gravity assist is an excellent technique to speed up or slow down a 
spacecraft, but is also a mechanism to deviate its course. For these reasons, NASA and ESA 
use the gravity assist as a helpful tool to reduce mission cost and to save fuel and travel 
time. Otherwise, planet exploration would be hard to pursue using current technologies. 

This phenomenon, also called slingshot effect or swing-by maneuver, can be easily under- 
stood as an elastic collision of two particles.-^ However, there is a paradox in the elastic 
derivation, because if one assumes that the planet's mass is much greater than the space- 
craft's, one would conclude that the initial and final speeds of the spacecraft (both when it 
is very far from the planet) should be equal, like in the usual Kepler motion. The "mistake" 
with this thinking is that one is forgetting the reference frame, since the above conclusion is 
correct for an observer on the planet, but not for an observer on the Sun due to the relative 
motion between the planet and the Sun. Therefore, one needs to make a Galilean addition of 
velocities for the latter case that results in the swing-by maneuver. One can check that the 
elastic collision is satisfied, the energy gained by the spacecraft is equal to the energy lost 
by the planet, but this is so small that the planet's motion is not altered.- The reader can 
also look for other approaches to explain the slingshot effect, like derivations from special 
relativity,- 1 ^ classical Lagrangian 5 and work- impulse.- Our goal in this paper is to get some 
insights of how the gravity assist changes the course of a spacecraft away from the ecliptic 
plane like in Ulysses mission. Our motivation is that we can find no literature about 3D 
gravity assist at the introductory physics level. 

Let us start with the description of our assumptions. We will use the word planet for the 
celestial body that alters the spacecraft's trajectory. 

1. For simplicity, the Sun, the Earth and the planet are all in the same plane, the ecliptic 
plane. We also consider that the Sun's and the planet's equators lie in that plane. 

2. The spacecraft motion can be divided in four parts, according to which celestial object 
contributes more to the gravitational force: 

(a) escaping the Earth, 

(b) moving by Sun interaction, 

(c) "swinging by" the planet and 
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(d) orbiting around the Sun or escaping the Solar System. 

We concentrate our work on the last two parts, where the gravity assist effect is 
employed and causes the final trajectory. 

3. After escaping the Earth, the spacecraft moves in a plane that is parallel to ecliptic. 

4. The interaction between the planet and the spacecraft is considered as scattering 
phenomenon in the planet's frame, occurring in a very short time compared to the 
period of the planet's circular orbit. Therefore, the planet's motion is taken as a 
straight line during the slingshot effect. 



II. HYPERBOLIC TRAJECTORY REVIEW 



Our goal in this section is to determine the scattering angle of the spacecraft, /3, due to 
the gravitational interaction with the planet, in the planet's frame. A hyperbolic trajectory 
is shown in Fig. HJ In polar coordinates, the equation for this kind of conic section is given 
by 

r^-C- 1 ', (1) 

1 +ECOS0 ' V ' 

where e(> 1) is the eccentricity and a is the semi-major axis. 




FIG. 1: Hyperbolic trajectory. F is the focal point where the force center is located. 
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As seen from this equation, a particular hyperbola is completely defined by two param- 
eters: e and a. In celestial mechanics, these parameters are usually written as functions of 
the first integrals of the motion: the energy, E, and the angular momentum, /, of a particle 
moving under the influence of a central inverse-square law force. In scattering problems, the 
impact parameter is used in place of the angular momentum. 

We consider that the initial velocity Ui and the pericenter r m are our known values, so 
that they will take the place of the above parameters, describing the spacecraft's hyperbolic 
orbit. 

By solving for the eccentricity in Eq. (CQ) when 6 = and r = r m , we find 

e = l + ^, (2) 
a 

where the semi-major axis is obtained by the usual celestial-mechanics equation 

GMm GM 

since E = \muf, when the spacecraft is far away from central force source. 

By taking r — > oo in the conic section equation (OQ) and solving for 9, we determine the 
angle between the asymptote and the x-axis 

cos# s = — - . (4) 

e 

With the geometry shown in Fig. [TJ the spacecraft's scattering angle is 

= 20s ~ 7T, (5) 

which has values < (3 < tt (attractive scattering). 



III. THE SIMPLE 3D CASE FOR GRAVITY ASSIST 

Let us first study the case when the spacecraft's- hyperbolic-orbit plane (SHOP) is per- 
pendicular to the ecliptic plane (Fig. [2]). Note that under assumption HI SHOP is a fixed 
plane in the planet's frame, but it is a constant-velocity moving plane in the Sun's frame. 
For this particular case, the 3D vector analysis gets simplified and allows us to picture how 
the general case can be established. We also choose that the spacecraft flies over the planet's 
south pole. Our goal here is to determine the angle between the spacecraft 's-final-orbit plane 
(SFOP) and the ecliptic plane, we will call it the elevation angle, 7. 
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FIG. 2: The simple 3D case for gravity assist: (a) its visual representation in 2D and (b) the top 
view. 

In a planet's frame as that shown in Fig. El the components of the planet's velocity 
relative to the Sun, Vp, are 

Vp,x = Vpcosa, Vp tV = 0, Vp >z = Vp sin a, (6) 

where a is the angle between the planet's velocity and the intersection of SHOP and the 
ecliptic plane. We consider < a < it, meaning that the spacecraft is approaching a planet 
outside the Earth's orbit. 

The initial and final velocities of the spacecraft in the Sun's frame, Vi and Vf, are obtained 
by Galilean addition of velocities 

Vi = Ui + Vp, v f = u f + V P , (7) 

where Vp obviously has the role of the velocity of the planet's frame with respect to the 
Sun's frame and Ui and itf are the spacecraft's initial and final velocities in the planet's 
frame. 

With the aid of Fig. [31 we find the magnitude V; L to be 

v\ = ( Ui + V P<X ) 2 + V£„ + V^ z = u 2 i+V^ + 2 Ui Vp cos a. (8) 
Note that a is also the angle between Vp and Ui. 
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FIG. 3: The spacecraft's-hyperbolic-orbit plane (SHOP) and the spacecraft's motion as seen from 
the planet. The vector Vp lies in the zx-plane. 

For the final velocity Vf, we first compute the components of Uf in the planet's frame 

Uf, x = Uf cos/5, Uf^y = Uf sin/3. (9) 

Therefore, the magnitude vj is given by 

v) = {u f , x + V P , X ) 2 + (u fi y + Vp, y f + V^ z = u 2 f + Vp l + 2u f V P cosa cos/3. (10) 

Then, we subtract Eq. ([H]) from Eq. ( jTUl) and obtain 

v j — v 2 = 2uVp cos a (cos/3 — 1), (11) 

since Ui = Uf = u by assumption HI 

Eq. ( TTTT) allows us to distinguish three situations according to how initial and final veloc- 
ities in the Sun's frame are related, 

• Vf > Vi\ The spacecraft increases its speed, if the planet scatters it {(3 ^ 0) and they 
encounter each other (| < a < n), 

• Vf < Vi\ It slows down, if the planet scatters it ((3 ^ 0) and the spacecraft "tries" to 
catch the planet (0 < a < |), 

• Vf = vf There is no speed change, if there is no scattering (/? = 0) or the spacecraft 
moves exactly perpendicular to the planet's motion (a = |). 
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FIG. 4: The spacecraft's- final-orbit plane (SFOP) in the Sun's frame. Vf t ± is like the lever arm. 

After the slingshot effect has occurred, we can determine the spacecraft's final orbit. 
According to assumption [2j the initial kinematic variables for the spacecraft orbiting around 
the Sun [part (d)] are those at the end of the slingshot [part (c)]. Thus, the initial velocity 
in part (d) is equal to Vf and the initial vector position r is taken to have a magnitude 
equal to the planet's distance from the Sun and lying in the ecliptic plane (Fig. [4]). This 
approximation can be justified because the gravitational force due to the planet is much 
weaker than that due to the Sun (see a numerical calculation for Jupiter at the end of 
Sec. ED . 

From the angular momentum conservation and the properties of the cross product, the 
SFOP is completely defined by Vf and r. Therefore, we conclude through Fig. 0] that 

tail7 = ^, (12) 

V f,z> 

where these velocity components are defined in the Sun's frame that has the x-axis aligned 
in the direction of the vector r. 

We recall Eq. (J7|) and use Fig. [5] to find 

tan 7 = 77 — = 777 , x n- ( 13 ) 

Vp + u f^ x cos ol (Vp/u) + cos a cosp 
Note that we have derived an expression that depends only on information before the sling- 
shot. 




FIG. 5: Top view for the velocity components of the spacecraft in the Sun's frame. 



IV. THE GENERAL 3D CASE FOR GRAVITY ASSIST 

Let us now study the case when the angle between SHOP and the ecliptic plane is some 
value 5, not necessarily 90°. We repeat the same derivation procedure done in the above 
section, but including now the effect of the angle 5. 




FIG. 6: The general slingshot case: (a) its visual representation in 2D and (b) the front view. 



From Fig. [61 the components of Vp in the planet's frame are 

Vp )X = Vp cos a, Vp tV = Vp sin a cos 5, Vp >z = Vp sin a sin 5, (14) 
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where a is still the angle between the planet's velocity and the intersection of SHOP and 
the ecliptic plane. 

Using Fig. [3] and Eq. ([7]), we find the magnitude as we did before 

v ] = u \ + Vp + 2uiV P cos a. (15) 

— * 

Note that a still represents the angle between Vp and Ui and the velocity components of the 
spacecraft in the planet's frame, U{ and Uf, do not change due to the SHOP orientation. 
The magnitude Vf follows from Eqs. (ED) and in Eq. (J7J) 

v j = uj + Vp + 2ufVp(cos a cos (3 + sin a sin (3 cos 5). (16) 

Then, we can subtract Eq. (I15p from Eq. fll6p and obtain 

vl — v\ = 2uVp (sin a sin (3 cos 5 + cos a (cos (3 — 1)) , (17) 

since Ui — Uf — u. With this result, the reader can study the ranges of a, (3 and 5 that make 
the spacecraft increase or decrease its speed or even those that lead to no speed change. She 
can also verify the usual 2D case (5 = 0, 7r).-^ 




FIG. 7: The velocity components in the Sun's frame: (a) front view and (b) top view. 



From the Fig. [7] and Eq. (fl2l) . the elevation angle is given by 



tan 7 — 

Vf, z > Vp + Uf tX cos a — u f, y cos o cos(7r/2 + a) 

sin (3 sin 5 



(Vp/u) + cos a cos (3 + sin a sin /3 cos 5 
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V. NUMERICAL CALCULATIONS 



Let us compute the elevation angle using Ulysses spacecraft data. From any college 
physics textbook,- we have the following Jupiter information: 



Mj = 1.90 x 10 27 kg, aj = 7.78 x 10 11 m, Rj = 6.99 x 10 7 m, tj = 3.74 x 10 8 s, 



where they are its mass, its distance from the Sun, its radius and its period, respectively. 
From Ulysses web page,- we have the spacecraft information: 



m = 366.7 kg, u = U{ = 13.896 km/s, Wj = 16.184 km/s, r m = 6.3 Rj. 
Therefore, the equations of Sec. HT1 result in (G = 6.67 x 10 -11 N • m 2 /kg 2 ) 



Assuming uniform circular motion, we estimate Jupiter's speed, Vp, and use Eq. (fT5|) to get 

a, 



We let the angle 5 be a free parameter. Table [J shows the final speed and the elevation 
angle computed from Eqs. ( fT6l) and (fl8l) with the above data. The usual condition for 



considered at the end of Sec. IIIII for deriving 7 (Ms = 1.99 x 10 30 kg). 

From Table [J we conclude that the minimum elevation angle that would keep the space- 
craft orbiting around the Sun is around 48°, meaning that the spacecraft would fly over 
a Jupiter's pole. Lower values of 7 would imply that the spacecraft will escape the Solar 
System. Notice that 5 represents the maximum latitude (both north and south) that SHOP 
crosses in the planet, under the assumption [1] (for 5 > 90°, the related parallel is 180° — 5). 
By symmetry, one could infer the results for 180° < 5 < 360°. When the spacecraft en- 
counters the planet and reaches the latitude 20.3° S as maximum (no further south), its 
final orbit around the Sun will be exactly perpendicular to the ecliptic plane. This situation 
repeats when it "tries" to catch the planet and reaches the latitude 20.3° N as maximum (no 
further north). The maximum speed given by gravity assist is obtained when the spacecraft 
"chases" the planet on the ecliptic plane and the minimum when they encounter each other 
on the ecliptic plane. 



a = 6.56 x 10 s m, e = 1.67, S = 127°, /3 = 74°. 




a,j = 13.1 km/s, a = 106°. 



elliptic orbit E < can be rewritten as v < ^/2GM s /r = 18.5 km/s, with r = aj as it was 



10 



5 (deg) 


Vf (km/s) 


7 (deg) 





26.0 


0.0 


15 


25.7 


8.0 


30 


25.1 


16.1 


45 


24.0 


24.1 


60 


22.5 


32.1 


90 


18.4 


48.0 


120 


13.0 


64.1 


146.9 


7.4 


80.0 


150 


6.8 


82.1 


159.7 


4.6 


90.0 


165 


3.5 


95.9 


170 


2.4 


104.5 


175 


1.4 


122.7 


180 


0.8 


180.0 



TABLE I: The final speed and the elevation angle as functions of 5. 

Let us check the agreement between our results and the actual Ulysses's orbit data. The 
Ulysses's elevation angle is around 80°, then from Tabled v = 7.4 km/s, which implies that 
its semi- major axis is about 3.10 AU (we have used Eq. ([3]) with E/m = v 2 /2 — GMs/aj). 
This conclusion is quite close with the actual Ulysses's semi-major axis 3.37 AU (computed 
by using the first aphelion 5.40 AU and the first perihelion 1.34 AU, obtained from the web 
page).- This comparison shows that our 3D slingshot formulation is numerically acceptable 
in spite of our assumptions and approximations. Moreover, it provides physical insight into 
how the gravity assist is used to deviate the spacecraft not only on the ecliptic plane but 
also away from it, according to the navigator's desire. 

Finally, let us numerically justify our approximation in Sec. Ill II for Jupiter case. The 
Sun's gravitational force around Jupiter is stronger than that of Jupiter when objects are 
located at distances of r > -J Mj / Ms aj = 0.03 aj = 334 Rj, so that our claim in Sec. II III is 
fairly good. 



11 



Acknowledgments 

I wish to thank my professor Daniel Azofeifa for getting me involved in the slingshot effect 
when I was a graduate student. I should also thank him for all discussions and comments, 
now that I have finally written down my thoughts. 



* Electronic address: hmorales@fisica.ucr.ac.cr 

1 A. A. Bartlett, and C. W. Hord, "The slingshot effect: explanation and analogies," Phys. Teach. 
23, 466-473 (1985). 

2 J. A. Van Allen, "Gravitational assist in celestial mechanics - a tutorial," Am. J. Phys. 71 (5), 
448-451 (2003), and references therein. 

3 W. Greiner, Classical Mechanics: point particles and relativity (Springer- Verlag, New York, NY, 
2004). 

4 J. J. Dykla, R. Cacioppo, and A. Gangopadhyaya, "Gravitational slingshot," Am. J. Phys. 72 
(5), 619-621 (2004), and references therein. 

5 K. J. Epstein, "Shortcut to the Slingshot Effect," Am. J. Phys. 73 (4), 362 (2005). 

6 C. L. Cook, "Comment on 'Gravitational slingshot'," Am. J. Phys. 73 (4), 363 (2005). 

7 R. Cacioppo, J. J. Dykla, and A. Gangopadhyaya, "Reply to 'Comment on Gravitational sling- 
shot'," Am. J. Phys. 73 (4), 363-364 (2005). 



Ulysses web pages: < jhttp: //ulysses . jpl .nasa. gov/| >, < |http: //u lysses -ops . jpl . esa. int/[ > 
The initial speeds Ui and vi are chosen from the ESA- JPL file: weekbody.txt, when Ulysses 
had its minimal speed in Jupiter's frame (date: December 16, 1991), as the frontier separating 
parts (b) and (c) of assumption [2j Strictly, m is the total mass at launch. 
9 See, for example, R. A. Serway and J. W. Jewett, Physics for scientists and engineers, 7th ed. 
(Thomson-Brooks/Cole, Belmont, CA, 2008). 



12 



